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OPERATOR  VALUED  FUNCTIONS  AND  BOUNDARY  VALUE  PROBLEMS 
FOR  THE  HELMHOLTZ  EQUATION.  I.  SPHERICAL  GEOMETRIES. 

1 . Introduction. 

Let  ft  be  an  open,  connected  region  in 
closed,  bounded  boundary  3ft.  We  denote  by 
where  ft  = ft  u 3ft. 

A number  of  boundary  value  problems  of  mathematical  physics 
reduce  to  requiring  solutions  of  the  Helmholtz  equation. 

LW  :=  (V2+k2)w  = f,  k2  > 0 (1.1) 

defined  in  the  region  ft  and  which  are  required  to  satisfy  certain 
boundary  conditions  denoted  typically  by  (be) . 

It  can  be  shown  that  the  problem  of  solving  the  boundary  value 
problem 

LW  = f,  w c (be)  (1.2) 

can  be  reduced  to  that  of  solving  a boundary  integral  equation  of 
the  form 

(I-K)u  = g (1.3) 

where  K is  a linear  integral  operator  mapping  H(3ft)  = H,  a 
Hilbert  space  of  functions  defined  on  3ft,  into  itself,  I denotes 
the  identity  operator  on  H whilst  u and  g are  respectively 
unknown  and  known  elements  of  H. 


E3  that  has  a smooth 
ft  the  closure  of 


-2- 


Strictly,  the  operator  K in  (1.3)  depends  on  the  parameter 
k appearing  in  (1.1)  and  we  should  in  consequence  regard  K as 
ar.  operator  valued  function  of  k.  It  will  be  our  intention  to 
make  explicit  this  aspect  of  the  boundary  integral  operator  K. 
Specifically,  we  consider  here  the  equation 

XKu  = u,  X e <t  (1.4) 

and  in  this  and  subsequent  communications  we  shall  concern  ourselves 
with  the  following  problems. 


(i)  The  determination  of  o(K),  the  spectrum  of  K. 

(ii)  The  determination  of  the  variation  of  o (K) 
with  k. 

(iii)  A discussion  of  the  significance  of  the  values  X = ±1. 

2 . A particular  case;  spherical  geometry. 

The  exterior  Neumann  problem  for  the  Helmholtz  equation  can  be 
reduced  to  a study  of  the  boundary  integral  equation  [Kleinman  and 
Roach,  1974] 

(I-K)w  = g (2.1) 

Here  K : I^Oft)  -*•  I^Oft)  is  a linear  operator  defined  by 

= f In- ' 

J 3ft  p q 


(Ku) (p) 


p,q  c 3ft 


(2.2) 
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where 


Y(P,Q)  = -elkR/2TrR,  R = R (P , Q)  , P,Q  e ft, 


i = n • V 

3np  p p 

A 

and  np  is  the  inward  drawn  unit  vector  normal  to  3ft  erected 
at  p e 3ft  . 

We  consider  the  equation  (2.1)  in  the  particular  case  when  ft  is 
the  unbounded  region  exterior  to  a sphere  of  radius  a.  Two 
situations  present  themselves;  k = 0 and  k ^ 0,  we  treat  each 
in  turn. 

2(a).  k = 0,  expression  for  eigenvalues. 

In  this  case  the  full  form  of  (2.1)  is 


“lp>  - d J,„  •<«>  I H-(l) 


- dsq  = g(p) , R = R(p,q) 


(2.3) 


The  expansion  of  1/R  in  terms  of  complete  sets  or  function  is  well 
known  (c.f.  Magnus  and  Oberhettinger , 1949/  which  hereinafter  we 
refer  to  as  M.O) . Specifically 


oo  r 


“ ■ Jo  ^ • V-  °> 


(2.4) 


P ;=  P(rp,0p,<j)p)  , 

Q :=  Q(rQ,6Q,<l>Q) 

P s=  P(rp' W ' 

q s-  q<rq,eq,*q) 

where 
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cos0  :=  cos  0:=  cos  0 cos  0 + sin  0 sin  0 cos(<b  -*  ) (2.5) 

pq  p q p q Yp  vq  ' 


r<  : = min(rp,r^) , r>  :=  max(rp,r^) 


(2.6) 


Consequently 


3 1 


3np|R 


n-1 

00  nr 

1 


n=0  r 


• P (cos  0) , 
n+1  n 


- I (n+1)  — 'n+2  pn(cos  0)  ' 


n=0 


r < r 
P q 


r > r 

p q 


(2.7) 


When  rp  = = a [i.e.  both  points  on  the  surface  9(2]  and  allowing 

for  approaching  from  D_,  the  inside  or  D+  the  outside,  of  9(2,  we 
take  the  mean  of  the  two  results  in  (2.7)  and  obtain 


3 1 


3np|R 


4 £ 4 • p„(cos 

n=u  a 


to  be  interpreted  as  a distribution  on  b (9(2).  The  homogeneous  form 
of  (2.3),  in  full  form,  can  now  be  written. 


W(p)  - ^ 


n />  2 fr  » 
d0  d<f>  a sin  0 

0 q io  q c 


(-  l P„(cOS  0)1  w(q)  = 0 

( 2a^  n=0  n > 


or 


\ 00  fit  f2n 

W(P)  + 47  E d6q  I d<P 
n=0  J0  q Jq 


sin  0 P (cos  0)w(q)  = 0 
q q n ^ 


(2.8) 


Solutions  of  this  equation  can  be  found  in  terms  of  the  spherical 


harmonics : 
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Y pn  (cos  Q)cos  m<]>  , e =»>  even  in  <p 

II  / k£  il 

Y™  n(0,4>)  :=  Pm(cos  0 ) sin  m<J>  , 0 =£>  odd  in  4> 

n,  u n 


where  P™(cos  0)  is  the  n-th  associated  Legendre  function  of  order 
m.  The  system  |y™  e,Y™  q|  is  known  to  be  complete  for  n 0, 
m £ n on  ( (0 , it)  x ( 0 , 2 tt  ) ) = L2(3ft)  [c.f.  M.O.  and  Smirnov,  Vol. 

IV,  1964].  If  in  (2.8)  we  set: 


w(p>  - *n,.‘VV  • 
0 

and  similarly  for  w(q)  then  (2.8)  becomes: 


(2.9) 


Ym  (0  ,<)>)+  t 
n, e p 4 


00  7T  r 2u 

n l d* 

J 1=0  >0  q Jo 


sin  0 P (cos  0 ) Y (0  , 4>  ) = 0 
q q n,e  q'^q' 

0 


This  double  integral  can  be  evaluated  by  the  orthogonality  proper- 
ties of  the  spherical  harmonics  [M.O.  p.  55]  and  we  obtain 


„m  i . , i A 4 ir  ..m  . . . 

n,e  0p'^p  47  2n+lYn,e^  p'^p 


Ym  (0 
n,e  p 


'V{1  + 5H+t}  ■ 0 


Consequently  from  2.10  we  see  that  the  eigenvalues  of  (2.3)  are  the 


scalars 


An  = - (2n+l) 


(2.11) 


and  *8  an  e^9envalue  of  (2.3)  of  multiplicity  2n+l  when 

n 0*  This  follows  because  for  fixed  n we  have  n eigenfunctions 
® — m — n for  the  even  situation  and  1 £ m £ n for  the  odd; 
hence  multiplicity  2n+l . 
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> 


n=2 

mult=5 


n=l 

mult-3 


n=0 

mult=l 


"0 1 &- 


-4 


-2 


4- 


-5 


-3 


-1 


A-plane 


Eigenvalue  distribution;  Afi(k) 

k=0 


2.  (b)  k / Q,  expression  for  eigenvalues. 

In  the  case  k ^ 0 the  homogeneous  form  of  equation  (2.1) 
becomes 

ikRi 


w(p)- 


IFj 


w (q) 

3ft  p 


dS  = 0. 

q 


(2.12) 


Using  known  expansions  (M.O.  p.  21]  we  have  for  this  case 
ikR  °°  (1) 


R = ik  (2n+l) jn (kr<)hn  ' ,(kr>)pn(cos  0) 


(2.13) 


where 


zn(X)  Zn+1/2(X) 


(2.14) 


and  represents  either  j the  spherical  Bessel  function 

°r  ^n^^  t*le  sP^er-ical  Hankel  function  of  the  first  kind  while 
Z represents  the  cylindrical  counterparts. 


t 


Proceeding  as  before  we  differentiate  (2.13)  with  respect  to 
n^,  obtain  expressions  for  p tending  to  3ft  from  inside  and 
outside  3ft  then  add  and  finally  divide  by  2 to  obtain: 


f^= — )=  l (2n+l)(j  (ka)h(1)  ' (ka)+j  ' (ka)h(1)  (ka)}p  (cos  G) 

3n  R 2 Ln  I n n Jn  n I n 

p ' n=0  v ' 


(2.15) 


Substituting  (2.15)  into  (2.12)  we  obtain: 


X f2lr  O iV^°°  f / 1 \ • 

W(p)-  de  d<j>  w(q)azsin6  J(2n+lMj  (ka)hvx;  (ka) 

271  JO  qJ0  q q 2 n=0  l n n 


..,2  » 


+ JA(ka>hn1)(ka>K 


(cos  0)  = 0, 


If  we  now  interchange  the  order  of  integration  and  abbreviate  the 
integrand  in  an  obvious  way  then  we  obtain: 


w(p)  - f ^ — y (2n+l)(j  h^'+j'h^l  d0  d^1  sinQ  P (cos  O)w(q) 
4tt  n£0  lJn  n Jn  n JJ0  qjQ  q q n 


As  before  set 


w(p)  = Y®  e (6  <J>  ) 

n»o  p r p 


= 0.  (2.16) 


(2.17) 


and  employ  the  associated  orthogonality  properties  [M.O.  p.  55] 
to  obtain 


CS(W-  Ii-(ka>2<2"+1»{jnhn1),+JAhn1)}{jSTT'’C,a}  - 0. 


This  implies  that 
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(0  ,<p  )(l-A.  (ka)2 
P P { 1 


[j  h*1* 
n n 


+h 


(1) 


0 


and  consequently  we  deduce  that 


A 

n 


i (ka) 2 [ j 'h(1)+j  h (1) '] 
n n n n 


(2.18) 


which  is  an  eigenvalue  of  (2.12)  of  multiplicity  2n+l(n^0). 

Now  it  can  be  shown  [Abramowitz  and  Stegun,  1964]  that  the 
Wronskian  is  given  by 


W ( j , h (^ ) ) : = J h(1)  ' 
n n n n 


i 

2 * 
(ka)  Z 


Hence  (2.18)  can  be  rewritten  in  the  form: 

n i (ka)  2 ^i/(ka)i+2jn'hiJ1)  } 


or 


A = 
n 


-l+2i (ka) 2 j 'h(1) 
Jn  n 


(2.19) 


where  it  is  to  be  understood  that  j 'h^  = j'(ka)h^(ka) 
and  ' denotes  differention  with  respect  to  the  argument.  Here 
the  eigenvalues  A^  are  dependent  on  the  value  of  k. 

Consequently  we  would  expect  that  as  k -►  0 we  would  recover  the 
eigenvalue  spectrum  obtained  in  the  previous  section.  That  this 
is  indeed  the  case  can  be  seen  by  means  of  the  following  asymptotic 
analysis.  From  [M.O.  p.  22)  we  have  the  expansions: 


(z) 


i-n-1 


n 

l <-Dm 

m=0 


(n+l/2;m) 

(2iz)m 


jnU)  - 


/tt 

“7 


. fz)ny  (~l)m(z/2) 2m 
l7Jm£0  n> ! T (n+m+3/2) 


(2.20) 


(2.21) 


r 


- 1 


V 


where 


and 


(n+1/2 ,m) 


(-l)m  (n+1)  (-n) 

m m 

m! 


(v) 


m‘ 


r (v+m) 
r (v) 


consequently  (2.20)  can  be  written  in  asymptotic  form  (for  small 
z)  as 


h^1 (2)  = i-"-1  S 


1Z 

z 


(-Dn  — — 1//2,^)  [ 1+0  (z)  ] 


(2iz ) 


xz 


I ~rwi+l  ^n+1/2  »n)  U+0  (z)  ] 
2 z 


and  we  can  deduce  that 


h:u'w»^  u+o«2„ 

i2  z 


Similarly  we  obtain 

_ /if  fzln  1 


III  TW2>  + 


and 


^ i / _ \ _ /iT  . nz 
Jn  ' Jn+T  f (n+3/2 ) 


n-1 


+ 0 (zn+1) . 


Consequently 


(1 ),„,/"/?  z"’1 


(2.22) 


(2.23) 


(2.24 


(2.25) 


jn (z)hn  (z)  = \7n+T  * 


r (n+3/2) 


+ 0 (zn+1 } " ;^n+n+l,n?  l} 

} 


= [anzn“1+o (zn+1) ] |bnelz/zn+1 [1+0 (z) ] 


iz 


= a b 
n n 


[1+0  (z) ) 


(2.26) 


t 
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where  a and  b are  defined  in  an  obvious  manner, 
n n 


Setting 


and 


c ; _ a b = n>,Tr  (n+1/2  *n)  = 1 n/rT (n+1/2 ,n) 

n n n 2n+1r (n+3/2) i2n  1 22n+1T (n+3/2) 


z = ka  , 


(2.27) 


we  find  on  using  (2 . 2 6 ) ,( 2 . 27)  that  (2.19)  can  be  written  in 
the  form 

1 


A = 
n 


2 ( elz 

-l+2iz^|cn^T[l+0(z)] 


(2.28) 


A = 
n 


-l+2ic  elz+0(z) 
n 


(2.29) 


In  order  that  this  result  (2.29)  should  agree  with  (2.11)  we 
must  have 

2n 

2 ic  = v - . 

n 2n+l 


(2.30) 


That  this  is  indeed  the  case  can  be  seen  as  follows.  Prom  the 
definition  of  (n+1/2, m)  and  (2.27)  we  obtain 

, n/?(n+l/2.m)  , nr U/2)  [-D n (n^l) n (-n) n 


2ic 


n 2^"r (n+3/2)  22n  nlT (n+3/2) 


and  since 


(-l)n(n+l)n(-n)n  = r(2n+l)  = 22nn! 


we  find  that 


2.  = nr (1/2) r (n+1/2) 


n 


n r (1/2) T (n+3/2)  n+1/2 


as  required. 
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For  large  values  of  z it  follows  from  (2.20)  and  the  fact  that 
(n+1/2,0)  = 1,  that  for  the  spherical  Hankel  function  of  the 
first  kind, 

hiJ1)(z)  = i‘n_1  ^ (1+0  (1/z ) ) 

and 

i z 

h(1)(z)  = i-n  - — (1+0 (1/z) ) . 
n z 


Similarly,  for  the  spherical  Hankel  function  of  the  second  kind 
we  have  that 


(z) 


^n+l) 


-iz 

~ — (1+0 (1/z ) ) . 


Hence 


and 


j„(z) 


(1) 


n 


(z)+h 


(2) 

n 


1 

2z 


. -n-1  iz . . n+1  -iz 
i e +i  e 

, 


(1+0 (1/z) ) 


5„(Z)C'<Z> 


.-2n+l  2iz 
i e +i 


(1+0 (1/z) ) . 


Therefore,  from  (21.9)  and  the  definition  of  the  Wronskian  it 
follows  that 


X 

n 


l+2iz2 j h(1) • 
Jn  n 


= jl+  (-1) ne2iz-lj 


1+0 (1/z) 


Consequently, 
values  of  z 


the  asymptotic  form  of  the  eigenvalues  for  large 
is 

Xn  = (-l)ne-2lz [1+0  (1/z) ] . (2.31) 


The  results  of  numerical  calculations  of  these  eigenvalues 
and  their  reciprocals  are  plotted  in  Appendix  II.  The  first 
eleven  eigenvalues  (equ.  2.18)  are  shown  as  functions  of  ka 


12 


illustrating  how  the  asymptotic  behavior  (equs.  2.29  and  2.31 
for  low  and  high  frequencies  respectively)  is  approached. 

3 . On  the  polar  decomposition  of  K. 

The  results  of  the  previous  section  indicate  that  as  k 
varies  the  eigenvalues  of  K become  complex.  Furthermore,  as  k 
increases  the  modulus  of  the  eigenvalues  tends  to  unity.  This 
suggests  that  K should  be  decomposable  either  as  the  sum  or  pro- 
duct of  operators,  one  of  which  is  unitary,  in  order  to  be  able  to 
exploit  the  behaviour  of  the  eigenvalues  of  K for  large  k. 

One  way  of  effecting  such  a decomposition  is  afforded  in  the 
so-called  polar  decomposition  of  operators  [Naylor  and  Sell,  1971]. 
Specifically,  given  a linear,  normal  operator  K:  H-*H,  where  H 
is  a separable  Hilbert  space  then  there  exists  a unique  decompo- 
sition of  K in  the  form 

K = UR  = RU  (3.1) 

where  R:  H-*-H  is  a positive,  self-adjoint,  linear  operator  and 
U:  H+H  is  a unitary  operator. 


Now  the  boundary  integral  equations  in  which  we  are  interested 
have  the  typical  form 


. 

(I-AK)w  = g e H, 

A e £ . 

(3.2) 

Therefore,  if 

K is  normal  then  (3.2) 

can  be  rewritten  either 

in  the  form 

* 

| 

(U_1-AR)w  = g , 

g = Ug 

(3.3) 

or  as 

i 

(R_1-AU) w = g , 

g - Rg 

(3.4) 

L 
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*■ 

Consequently  in  the  light  of  the  results  of  the  previous 
sections  we  are  led  to  make  the  following  conjectures. 

Kl . As  z:  = ka-*-»  then  K->U  therefore  R-*±I  and  R ^-*-±1 

where  I is  the  identity  on  H.  Consequently,  the  high 
frequency  behavior  of  the  required  solutions  to  (3.2)  is 
determined  from  the  solutions  of  the  equation 

(I-XU)w  = g , w = lim  Rw  = ilim  w , g = lim  g (3.5) 

z~+°°  z-+°°  k -*■<*> 

where  this  equation  has  been  obtained  from  (3.2)  and  (3.4) 
in  the  limit  as 

K2 . As  z:  ka-*-0  , K becomes  self  adjoint  therefore  U-*-±I 
and  U ■'■->■+ 1.  Consequently,  the  low  frequency  behavior 
of  the  required  solution  to  (3.2)  is  determined  from  the 
solutions  of  the  equation 

( I-AR) w_  = g , w = lim  Uw  = lim  w , g = lim  g (3.6) 

° ° z+0  z+0  ° k+0 

where  the  equation  has  been  obtained  from  (3.4)  in  the 

limit  as  z+0. 

The  Spectral  Theorem  (Naylor  and  Sell,  1971)  states  that  if 
K:H+H  is  a compact,  normal  operator  then  there  exists  a resolution 
of  the  identity  anc*  a sequence  of  complex  numbers 

such  that 

K • I „nEn  (3.7) 

n 


where  the  convergence  in  (3.7)  is  understood  to  be  in  terms  of 
the  uniform  operator  norm  topology. 


r 


_! 
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{ } is  a resolution  of  the  identity  provided 

2 

(i)  are  orthogonal  projections,  i.e.,  En 

and  R(En)iN(En) 


E = E* 
n n 


EE  = 0 , n / m 
n in 


(ii) 

(iii)  l En  = 1 • 


In  the  case  under  consideration,  u and  E may  be 

n n 

determined  from  2.12,  2.15,  2.18  to  be 

v„  - A'1  = iz2{j_;(z)h|;ll(2)tjn(z)h|;i)'(2)} 

+1  f77 

77  Jo 


En(,)  = I <36 


q I 


2tt 

d4>  sin0  P (cosG  ) ( • ) . 
0 q q n p , q 


(3.8) 


(3.9) 


The  well-known  properties  of  the  spherical  harmonics  (M.O.  p.  55) 
enable  one  to  verify  that  (i)-(iii)  are  fulfilled. 

In  general,  provided  |yn|  ^ 0,  and  since  the  E^  are 
orthogonal  projections,  we  can  write 


K = l uE  = l |u  |E  I H Ip  I 1E  = RU 
“ n n n n n ^ m1  m'  m 


(3.10) 


Now,  any  operator  K which  can  be  written  in  the  form  (3.7)  is 
called  a weighted  sum  of  projections.  Therefore  the  operators 
R and  U appearing  in  (3.8)  and  defined  by 

-1. 


l KlEn  ' U:=  l ^nK1"  En 


(3.11) 


are  both  weighted  sums  of  projections  and  R is  self-adjoint 
since  the  weights  are  real  (|y  |eR).  Furthermore,  by  first 
recalling  the  properties  of  the  resolution  of  the  identity, 
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elementary  manipulation  shows  that  u is  unitary.  Therefore 
(3.10)  provides  the  required  polar  decomposition  of  K.  We  are 
now  in  a position  to  examine  the  conjectures  Kl  and  K2  if  we 
recall  the  asymptotic  estimates  obtained  in  the  previous  sections, 
namely 


A (z)  = < 

^-(2n+l)  + 0(z) 

, for 

small 

z 

and 

fixed  n. 

n 

[(-l)ne~2lz+  0 (l/z) 

, for 

small 

z 

and 

fixed  n. 

Consequently 

u (z) :=  A"1  (z) 
n n 


tiktit  11+0 <*>> 

(-l)ne2lz [1+0 (l/z) ] 


for  small  z and  fixed  n, 
for  large  z and  fixed  n, 

(3.12) 


and  we  obtain  by  direct  substitution  of  these  quantities  into 


(3.11 


J f P (cos  0 ) ( • ) dS  [ 1+0  ( z ) ] 

n=0  47ra2J9fin  pq  <1 


U = l 


-I (1+0 (l/z)  for  large  z . 

00  f 

- I Until  I P (cos  0)  (-)ds  (l+0(z)  ] 
n=0  47ra2  hti  n pq  q 

-l(l+0(z))  for  small  z . 


(3.13) 


(3.14) 


4 . Conclusion , 

In  this  note  we  have  begun  an  investigation  into  both  the 
theoretical  and  practical  consequences  of  treating  certain 
commonly  occuring  operators  of  mathematical  physics  as  operator 
valued  functions  of  a parameter. 

As  a first  consideration  of  this  problem  we  have  confined 
our  attention  to  boundary  value  problems  associated  with  the 
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4. 


Helmholtz  equation  in  the  particular  case  of  spherical  symmetry . 

We  have  seen  that  when  8Q  is  a sphere  then  equation  (1.4)  can 
be  solved  analytically.  Consequently  a number  of  quantities  which 
are  known  to  depend  upon  a certain  parameter  can  be  computed 
directly  and  so  the  dependence  can  be  exhibited.  Specifically, 
we  have  seen  that  as  k varies , the  eigenvalues  of  K become 
complex.  Furthermore  we  have  seen  that  as  k increases  the 
modulus  of  each  eigenvalue  tends  to  unity.  This  suggests  that 
K should  be  decomposable  as  the  product  of  operators,  one  of 
which  is  unitary,  in  order  to  exploit  the  behavior  of  the  eigen- 
values of  K for  large  k.  This  decomposition  was  explicitly 
exhibited  for  the  spherical  boundary.  Promising  as  this  appears 
at  first  sight,  it  turns  out  not  to  be  a very  practical  proposi- 
tion in  general,  since  its  implementation  requires  an  intimate 
knowledge  of  all  the  eigenvalues  and  eigenvectors  of  K.  Further- 
more such  polar  decomposition  of  K is  only  conveniently 
available  for  compact  normal  operators  and  while  it  can  be  shown 
(Appendix  I)  that  K is  normal  for  spherical  geometries,  this 
may  not  be  true  in  general.  Nevertheless  we  shall  show  in 
subsequent  communications  that  there  are  factorization  methods 
which  do  offer  good  practical  prospects  especially  if  instead  of 
K attempts  are  made  to  factorize  (I-K)  the  resolvent  of  K. 

Finally  we  remark  that  the  equation  whose  solution  is  sought, 
(2.1), has  obvious  difficulties  for  those  values  of  ka  for 
which  An(ka)=l,  where  A^(ka)  are  the  eigenvalues  of  equ.  1.4. 

To  show  that  such  critical  values  do  occur  in  the  spherical 
example  treated  here,  the  quantity  | (k.a)  —1 1 is  plotted  as 

functions  of  ka  for  0 < n < 10  in  Appendix  II. 


1 
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Appendix  I.  Concerning  the  normality  of  boundary  integral 
operators . 


On  L2(aft)  we  define  the  inner  product 

(U|V)  = u ( p ) v ( p ) dS  , p e a n 

) an  p 


(Al.l) 


where  dS  denotes  an  element  of  area  of  the  surface  3ft  taken 
P 


with  respect  to  the  coordinates  of  a point  p e a ft, 


We  introduce  an  integral  operator  K:  L2(9ft)-*  L2  (3ft) 


defined  by 


(Ku)  (p)  : = 


ay  (p,q) 


afiq  3np 


u (q) dSq  , p,q  e 3 ft 


(A1.2) 


where  the  kernel  of  this  operator  is  generated  by  Y(P,Q), 

P,  Q e ft,  a function  'of  position  in  ft  = ftU3ft.  The  domain 
functional  is  conveniently  taken  to  be  the  fundamental  solution  of 
the  associated  partial  differential  equation.  In  this  example, 
in  which  we  are  considering  the  Helmholtz  equation,  we  take 

_ ikR 


Y (P/Q)=  -• 


2ttR 


R:  = PQ  , P,Q  e ft  . 


(A1.3) 


The  adjoint  of  K is  obtained  by  considering 
(Ku,v)  = [ (Ku)  (p)  v (p)  dS 

hn  P 


= f •[  Ij-  (p,q)u(q)dS  v(p)dS 
hn  J an  3nn  q 


3fip;9qq  P 


q 

(u, K*v) 


3y 


= f u(q)‘  f ji—  (p,q) V (p)dS 
J»n  J9fip  3np  P 


dS. 
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where 


(K*v ) (q) : = 


Jn~  <p.q)v(P)ds  . 
3S1  p 
P 


(A1.4) 


In  order  to  investigate  the  normality  of  the  operator  K we 
examine  the  properties  of 

S:  = K*K  - KK* 

Substituting  for  K and  K*  in  Al-5  we  obtain 


(A1.5) 


(Su) (p) 


" I I u 

i SO  j ao 


(q) 


3Qq 


^~(r'q)lH“(r'p)"S_(p,r)^~(q,r) 


3rT 

r 


3n 


3n 


Consequently 

D (p,q)  : = 


|j-<r,q)|j-(r.P)-|^(P,r,|^(q,r) 


dS  dS 

q r 


(A1.6) 


dS  = 0 (A1.7) 

r 


is  a necessary  and  sufficient  condition  for  the  normality  of  K. 

To  see  that  Al.7  is  satisfied  for  the  Helmholtz  equation 
and  3f2  a sphere  we  first  write  (2.15)  in  the  form 


jj-'p.q)  - J0Pn(alPn(coS  0pq):-  JnFn(a)Pn(p,q) 


(Al.7) 


n=0 


ik2 


where  F (a)  = (2n+l)(j  (ka) h (1)  ’ (ka)  + j ' (ka)  h (1)  (ka)  ) 

n i Jn  n n n 

Substituting  Al.7  into  Al.6  yields 
/ [ 00  00 

D(P*q)  - l F (a) P (r ,q)  \ F~(a)Pm(r,p) 

' 3fir (n=0  n n m=0  m m 


- I F (a) P (p, r)  l F (a)P  (q,r) 
n=0  n n m=0  m m 


dS 
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“ I l Fn (a) F {pn(r,q)P  (r,p)-P  (p,r)P  (q,r)ldS 
n=0  m=0  n m Jafl  l n m n m M / r 


oo  • 

■ J0FnUlfn  (a)  1ST  {P„(p-q)'P„(P-qlJ 


= 0. 


Hence  for  the  particular  case  of  the  Helmholtz  equation  and 
a spherical  surface  the  boundary  integral  operator 
K:  1*2  Oft)  -*■  L2(3fi)  is  normal. 
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Appendix  II.  Numerical  results. 


Equation  (2.19)  defines  A^  explicitly  as 


X = ± 

n -l+2i(ka)2ji;(ka)hi51)  (ka) 


To  illustrate  the  behavior  of  the  eigenvalues  as  functions 

of  ka,  A^  is  plotted  as  a function  of  ka  for  0 n 10 

and  0 <_  ka  10  in  Figures  1-11.  Notice  that  the  asymptotic 

behavior  for  small  and  large  ka  is  evident  although  for  the 

larger  values  of  n,  values  of  ka  greater  than  10  are 

required  before  the  asymptotic  approximation  is  reasonable.  In 

Figures  1-6,  8,  10  and  11,  the  horizontal  and  vertical  scales 

are  different  so  extra  care  is  urged  in  interpreting  these  graphs, 

The  reciprocal  eigenvalues  p = y—  are  plotted  in  Figs.  12-22. 

n n 

The  advantage  in  plotting  in  this  way  is  that  the  wn  are 
bounded  and  the  same  scale  can  be  employed  for  all. 

Since  the  equation  to  be  solved  is 
(I-Ak)w  = f 

with  A = 1,  it  is  of  interest  to  observe  how  close  this  value 
of  1 is  to  an  eigenvalue.  To  exhibit  this  the  quantity 
|A  (ka)-l|  is  plotted  as  a function  of  ka  for  0 <_  n <_  10 
and  0 <_  ka  <_  10  in  Figures  23-33. 


* 
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Figure  3.  An(ka),  n = 2 
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X-plane 


Figure  8.  A (ka) 
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y-plane 


Figure  19.  y (ka)  = , n = 7 
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Figure  29.  |An(ka)-l|,  n = 6 
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Figure  30.  |Xn(ka)-l|,  n = 7 
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